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Abstract 

We performed ^'^C0{1 - 0), ^^C0{1 - 0), and HCN(1 - 0) single-dish observations (beam size ~ 
14" — 18") toward nearby starburst and non-starburst galaxies using the Nobeyama 45 m telescope. The 
""^■^00(1 - 0) and HCN(1 - 0) emissions were detected from all the seven starburst galaxies, with the 
intensities of both lines being similar (i.e., the ratios are around unity). On the other hand, for case 
of the non-starburst galaxies, the ^'^C0(1 - 0) emission was detected from all three galaxies, while the 
HCN(1 - 0) emission was weakly or not detected in past observations. This result indicates that the 
HCN/^'^CO intensity ratios arc significantly larger (~ 1.15 ±0.32) in the starburst galaxy samples than 
the non-starburst galaxy samples (< 0.31 ±0.14). The large-velocity-gradient model suggests that the 
molecular gas in the starburst galaxies have warmer and denser conditions than that in the non-starburst 
galaxies, and the photon-dominated-region model suggests that the denser molecular gas is irradiated by 
stronger interstellar radiation field in the starburst galaxies than that in the non-starburst galaxies. In 
addition, HCN/^'^CO in our sample galaxies exhibit strong correlations with the IRAS 25 /xm flux ratios. 
It is a well established fact that there exists a strong correlation between dense molecular gas and star 
formation activities, but our results suggest that molecular gas temperature is also an important parameter. 
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1. Introduction 

Galaxies consist of vast numbers of stars formed from 
molecular clouds. It has therefore been suggested that 
the process of star formation in galaxies depends largely 
on the distribution, kinematics, and physical conditions 
of the molecular gas. One of the primary molecular gas 
properties critical to star formation is the number den- 
sity. Observational studies of molecular clouds in our 
Galaxy suggest that stars arc formed from the dense lo- 
cale of molecular clouds, rather than the diffuse regions 
(e.g., Lada 1992). Extragalactic single-dish molecular gas 
observations indicate that the overall amount of dense 
molecular gas, which can be traced by the HCN(1 - 0) 
luminosity, is tightly correlated with the global star forma- 
tion rate, which can in turn be traced by the far-infrared 
(FIR) luminosity observed with the Infrared Astronomical 
Satellite (IRAS) (Solomon et al. 1992; Gao & Solomon 
2004a; Gao & Solomon 2004b). This correlation is much 
tighter than that between the diffuse molecular gas, or in 
other words, the total amount of molecular gas that can be 
traced by ^^C0(1 - 0), and the global star formation rate 



(e.g., Sanders et al. 1991; Young & Scoville 1991; Young 
et al. 1996; Gao & Solomon 2004b). In addition, Solomon 
et al. (1992) and Gao & Solomon (2004b) also raise the 
possibility of the fact that the HCN/^^CO luminosity ra- 
tio correlates well with the FIR/^^CO luminosity ratio, 
namely the fraction of dense molecular gas correlates well 
with the star formation efficiency, and this suggests that 
HCN/^^CO is a good indicator for star formation. 

High spatial resolution interferometric observations to- 
ward the center of a galaxy with a high infrared luminosity 
indicate that extragalactic star forming regions coincide 
well with the distributions of dense molecular gas (HCN) 
(Kohno et al. 1999). The distributions of diffuse molecular 
gas (^^CO), on the other hand, only loosely match those of 
star forming regions. As a matter of fact, the HCN/^^CO 
integrated intensity ratios are high (~ 0.1 — 0.2) in ac- 
tive star forming regions in the centers of galaxies with 
high infrared luminosity (> 10^° L©), but low (< 0.1) in 
the centers of galaxies with low star forming activities or 
post-starburst activities (Kohno et al. 1999; Reynaud & 
Downes 1999; Kohno et al. 2002). Star forming regions 
in the outer parts of galaxies, where normally have low 
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infrared luminosities or low star formation rate, have low 
HCN/^^CO at giant molecular cloud scale (Heifer & Blitz 
1997; Brouillet et al. 2005). These results strongly imply 
that the density of molecular gas play a critical role in the 
active star formation phenomena. 

The relation between star formation and temperature, 
another important molecular gas property, has so far not 
been well established. There are several surveys observ- 
ing 12C0(2 " 1) and ^^00(1 - 0) lines toward galaxies 
with various star formation activities (such as starburst, 
interacting, and quiescent galaxies), and most of the sur- 
veys yield an average ^^C0(2 - 1)/(1 - 0) of around unity 
(Braine et al. 1993; Aaho et al. 1995; Hafok & Stutzki 
2003). Surveys focusing on the ^^00(3 - 2) and ^^C0{1 
- 0) lines indicate that variations in ^^C0(3 - 2)/(l - 0) 
tend to have somewhat larger dependence on star forma- 
tion activities as compared to ^^C0(2 - 1)/(1 - 0); star- 
burst galaxies tend to have higher ^^C0(3 - 2)/(l ~ 0) of 
around unity, but normal galaxies tend to have lower ra- 
tio of around 0.5 or lower (Hafok & Stutzki 2003; Dumke 
et al. 2001; Devercux et al. 1994). These multi- transition 
^^CO observations suggest that the ^^CO lines in nearby 
normal and starburst galaxies are optically thick, with a 
tendency to be saturated to unity, namely thermalized, 
in the vicinity of certain star formation activities. It is 
therefore difficult to deduce correlations between star for- 
mation activities and multi-transition ^^CO ratios. 

The intensities of the ^^C0(1 - 0) lines, on the other 
hand, display more drastic changes than those of the ^^CO 
lines with respect to star formation activities. In merg- 
ing galaxies, the ^^C0(1 - 0) intensities tend to be very 
low, and 12C0(1 - 0)/i^^CO(l - 0) often exceed 20, in 
contrast to normal or starburst galaxies (not including 
extreme starbursts) that have an average value of ^ 10 
(Aalto et al. 1991; Casoh et al. 1992; Aalto et al. 1995). 
However, this clear contrast in ^^00(1 - 0)/i3CO(l - 
0) is only manifested between merging galaxies and nor- 
mal/starburst galaxies, and there exists only a weak trend 
between normal and starburst galaxies as a function of 
dust temperature or the 60 fim and 100 fim FIR flux ra- 
tio (Sage & IsbeU 1991; Aaho et al. 1991; Aaho et al. 
1995). 

In the study of M51 multi-line observations, we dis- 
covered an excellent density and temperature probe, the 
HCN(1 - 0)/"CO(1 - 0) ratio (Matsushita et al. 1998, 
hereafter Paper I; Matsushita et al. 1999). The utility of 
this ratio has been confirmed by higher-J interferometric 
^^CO observations (Matsushita et al. 2004). 

In this paper, we present HCN(1 - 0), ^^00(1 - 0), 
and ^^C0(1 - 0) observations toward galaxies with various 
levels of star formation activities, and discuss the relation 
between the physical conditions of molecular gas and the 
levels of star formation activities. In addition, we compare 
the IR properties of these galaxies, and show that the 
molecular gas and dust are localized in similar regions 
and strongly irradiated by the radiation from starburst 
regions. 
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Fig. 1. Correlation diagrams between the 

/(12 iim)/f{25 fim) and the /(60 Atm)//{100 fj-m) IR 
flux ratios. Filled and open circles denote the starburst and 
the non-starburst galaxy samples, respectively. Our sample 
galaxies extend over a wide range of IR flux ratios reported 
in Helou (1986) and Dale et al. (2001), covering a wide 
spectrum of star formation activities. 

2. Sample Galaxies 

For this survey, we picked 7 starburst galaxies and 3 
non-starburst galaxies so as to cover a wide range of 
star formation activities. We summarized the basic pa- 
rameters of the starburst and non-starburst galaxy sam- 
ples in Table 1 and Table 2, respectively. To ensure the 
wide spectrum of star formation activities of our samples, 
we plot the IR fiux ratios between 12 ^m and 25 /im, 
/(12 /im)//(25 /im), and that between 60 jim and 100 /im, 
/(60 /im)//(100 /im) in Fig. 1. It has well established that 
this correlation traces star formation activities (Helou 
1986; Dale et al. 2001). Our sample galaxies extend over 
a wide range of the IR flux ratios, covering most of the 
range of the diagrams shown in Helou (1986) and Dale et 
al. (2001), thus ensuring the wide coverage of star forma- 
tion activities. 

2.1. Starburst Galaxy Sample 

Among the starburst galaxy samples are M82, 
NGC 253, NGC 2146, NGC 3504, NGC 2903, NGC 6946, 
and NGC 6951, whose FIR luminosities are on the order 
of 10^° Lq (Sanders et al. 2003). 

M82 and NGC 253 exhibit energetic phenomena at 
various wavelengths, such as Ha, X-ray, and molecular 
gas outflows (e.g., Strickland et al. 2004; Nakai et al. 
1987) and molecular superbubbles (e.g., Matsushita et al. 
2005; Sakamoto et al. 2006). These galaxies are therefore 
known as prototypical starburst galaxies (e.g., Rieke et al. 
1980; Rieke et al. 1988; Devereux 1989). 

NGC 2146 and NGC 3504 display star formation ac- 
tivities comparable to M82 and NGC 253 on the basis of 
1.65, 2.2, and 10 /im photometric observations (Devereux 
1989; Devereux et al. 1994). In particular, NGC 2146 
hosts a kpc-scale starburst-drivcn X-ray outflow (Armus 
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Table 1. 


Basic parameters 


of the starburst g 


:alaxy sampl 


es. 








TvDe* 


Q;(B1950) 


(5fB1950l 




Dist.* 


Ref.§ 


—J p 

-^FIR 


NGC 9^"^ 


SARl'sV 


nn'M5™05'^fi3 

V7V7 4nJ \jO . UO 




\ 


3 ^ 


O 




NGC 2146 


SB(s)ab pec 


06'n0"'41?10 


+78°22'28."1 


2 


17.2 


9 


10.93 


NGC 2903 


SAB(rs)bc 


09'^29'"20?26 


+21°43'20."74 


3 


8.9 


10 


10.05 


M82 


10 sp 


09''51"43?5 


+69°55'00."0 


4 


3.9 


11 


10.61 


NGC 3504 


(R)SAB(s)ab 


lli^00™28?53 


+28°14'31"2 


5 


26.5 


9 


10.56 


NGC 6946 


SAB(rs)cd 


20'^33'"49?2 


+59°58'49"5 


6 


5.9 


12 


10.01 


NGC 6951 


SAB(rs)bc 


20'^36'"36?59 


+65°55'46."0 


7 


26.8 


13 


10.46 



Notes. 

* Morphological type taken from RC3. 

t References for the positions. 

1 Distance in Mpc. 

§ References for distance. 



I Far infrared luminosity taken from Sanders et al. (2003). 

References — (1) Keto et al. (1993). (2) Condon et al. (1982). (3) 2.2 /im peak of Wynn-WiUiams & Becklin 
(1985). (4) Shen & Lo (1995). (5) Condon et al. (1990). (6) Turner & Ho (1983). (7) Saikia et al. (1994). 
(8) Rekola et al. (2005). (9) TuUy (1988). (10) Drozdovsky & Karachentscv (2000). (11) Sakai & Madore 
(1999). (12) Karachentsev et al. (2000). (13) Valentini et al. (2003). 



et al. 1995; Delia Ceca et al. 1999; Inui et al. 2005) and 
molecular outflow and superbubbles (Tsai et al. 2009) 
from the central region, which are similar to those in M82 
or NGC 253, and therefore strongly support the existence 
of ongoing starburst. 

NGC 2903 is one of the Sersic-Pastoriza galaxies (a cat- 
alogue of galaxies with peculiar or complex nuclei such as 
"hot spots" or "amorphous nucleus"; Sersic & Pastoriza 
1965; Sersic & Pastoriza 1967; Sersic 1973) and its nuclear 
appearance is classified as "hot spots" . This nuclear "hot 
spots" region is a powerful IR source, and it has been 
proposed that this phenomenon could be the result of a 
starburst (Rieke et al. 1980; Telesco & Harper 1980), a 
conjecture that has been supported by radio and IR obser- 
vations (Wynn-WiUiams & Becklin 1985; Alonso-Herrero 
et al. 2001). 

NGC 6946 is a large spiral galaxy with a small nuclear 
bar (Elmegreen et al. 1998) and ongoing nuclear starburst 
(Turner & Ho 1983; Engelbracht et al. 1996). 

NGC 6951 is also a Sersic-Pastoriza galaxy, with a cir- 
cumnuclear "hot spots" ring (Barthet al. 1995; Elmegreen 
et al. 1999). This ring is bright in Ha (Wozniak et al. 
1995; Gonzalez-Delgado et al. 1997) and exhibits high star 
formation rate (SFR) of ^ 3 — 4 Mq yr~^ as well as high 
star formation efficiency (SEE), comparable to those in 
the central regions of nearby starburst galaxies (Kohno et 
al. 1999; Wozniak et al. 1995). 

2.2. N on- Starburst Galaxy Sample 

The non-starburst galaxy sample consists of NGC 4736, 
NGC 4826, and NGC 5195, whose FIR luminosities are 
around 10^ (Sanders et al. 2003). 

The optical spectra of the nuclear regions of both 
NGC 4736 and NGC 5195 exhibit strong Balmer absorp- 
tion lines (Ho et al. 1995; Taniguchi et al. 1996; Sauvage 
et al. 1996; Greenawalt et al. 1998). Since these strong 
Balmer absorption lines are typically found in the spectra 
of A-type stars, stellar populations of galaxies with such 
absorption lines are believed to be dominated by A-type 



stars. Indeed, the stellar population synthesis analysis im- 
plies that the dominant stellar populations of NGC 4736 
and NGC 5195 extend up to at most A4 to A7, and that 
the optical light is dominated by these stellar populations 
(Pritchet 1977; Warner 1974). The likely mechanism for 
the creation of "A-type star dominated" galaxies is as fol- 
lows: Starburst phenomena produce enormous numbers 
of early to late-type stars, assuming a general initial mass 
function. Early-type stars (OB stars) live at most 10® 
years, so that at the time 10^ years after the beginning 
of the starburst, the early-type stars are burned out, with 
vast numbers of A-type stars remaining. Therefore, galax- 
ies with numerous A-type stars can be considered post- 
starburst galaxies (Rieke 1988; Rieke et al. 1988; Walker et 
al. 1988). Molecular gas observations support the hypoth- 
esis that, in the Balmer absorption dominated regions, 
there is little molecular gas (Tosaki & Shioya 1997). In 
some cases, molecular gas exists in the "A-type star domi- 
nated" galaxies, but it is stable against gravitational insta- 
bility (Shioya et al. 1998; Kohno et al. 2002) with almost 
no high-density gas present (Kohno 1998; Kohno et al. 
2002). Therefore, we classified NGC 4736 and NGC 5195 
as "post-starburst" as in other articles. 

NGC 4826 is classified as a type 2 transition object (the 
[O l] line from the nucleus has a strength intermediate 
between those of Hll nuclei and LINERs; Ho et al. 1997). 
Even if all of the Ha emission of this galaxy originates 
in the star forming regions, the extinction corrected total 
Ha luminosity, i(Ha), of 5.7 x 10^*^° erg s'^ (Young et al. 
1996) indicates a SFR of only ~ 0.45 Mq yr'^, using the 
relationship between the SFR and the L(Ha) (Kennicutt 
1998). Hence this galaxy can be safely regarded as a non- 
starburst galaxy. 

3. Observations 

Observations of the central regions of the sample galax- 
ies using the Nobeyama 45 m telescope were carried out in 
the HCN(1 - 0), 13C0(1 - 0), and 1^00(1 - 0) lines (rest 
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Table 2. Basic parameters of the non-starburst galaxy samples. 



Galaxy 


Type* 


a(B1950) 


(5(B1950) 


Ref.t 


Dist.t 


Ref.§ 


Activity II 




NGC 4736 


(R)SA(r)ab 


12M8™31?910 


+41°23'31"78 


1 


4.3 


3 


Post-SB 


9.59 


NGC 4826 


(R)SA(rs)ab 


12i>35™16?07 


+21°57'13"5 


2 


4.1 


3 


T2 


8.98 


NGC 5195 


10 pec 


13^27™53?27 


+47°31'25"5 


2 


8.4 


4 


Post-SB 


9.44: 



Notes. 

* Morphologieal type taken from RC3. 

t References for the positions. 

t Distance in Mpc. 

§ References for distance. 



" Star formation activity. "SB" and "T2" denote starburst and type 2 transition object, respectively. See text for details. 
* Far infrared luminosity taken from Sanders et al. (2003). A value with a colon (":") indicates that the value has a large 
uncertainty. 

References — (1) Turner & Ho (1994). (2) Hummel et al. (1987). (3) Tully (1988). (4) Feldmeier et al. (1997). 



frequency = 86.632 GHz, 110.201 GHz, and 115.271 GHz, 
respectively). We observed ^^C0{1 - 0) and HCN(1 - 
0) simultaneously for the starburst samples during 1997 
February - April, 1997 December, and 1999 March - April 
(the 1999 observation run was only for NGC 6951). The 
^^C0(1 - 0) and ^^C0(1 - 0) simultaneous observations 
for the non-starburst samples and M82 were performed 
during 1998 March - April. 

Two SIS receivers, SlOO and S80, were used, and these 
were tuned to the ^^C0{1 - 0) and HCN(1 - 0) lines, 
respectively, for the ^^CO and HCN simultaneous ob- 
servations. The half power beam widths (HPBW) for 
the ^•^CO and HCN simultaneous observations (except for 
NGC 6951) were 15" and 17", and the main-beam efficien- 
cies, ?7mb, were 0.51 ±0.03 and 0.47±0.02, respectively. In 
case of NGC 253, we applied T^mb of 0.48 and 0.35, respec- 
tively, due to the low declination (and therefore the low 
elevation at the telescope site) of this source. The HPBW 
for the NGC 6951 observations were 15" and 18", and 77,„b 
were 0.48 ±0.02 and 0.49 ±0.02, respectively. The differ- 
ences in the HPBW and rj^b from other observations are 
due to the different observation year for this galaxy. 

For the ^^C0{1 - 0) and ^^C0{1 - 0) simuhaneous ob- 
servations, SlOO and S80 were tuned to ^^CO and "CO, 
respectively. Since the different (S80 instead of SlOO) re- 
ceiver was used for the ^^CO line, the HPBW and Tymb 
were different from the ^'^CO and HCN simultaneous ob- 
servations; the HPBW for the "CO and ^^CO simul- 
taneous observations were both 14", and the Tymb were 
0.40 ± 0.03 and 0.51 ± 0.03, respectively. 

The telescope pointing was checked and corrected ev- 
ery hour, and the absolute pointing errors were less than 
4" during the observations. Due to the use of two SIS re- 
ceivers, there were no relative pointing offsets between two 
simultaneously observed lines. The single side-band sys- 
tem noise temperatures at 87 GHz, 110 GHz, and 115 GHz 
throughout the observations were 200 - 400 K, 300 - 
1400 K, and 600 - 1800 K, respectively. As back-ends, we 
used 2048-channel wide band acousto-optical spectrom- 
eters (AOS), with a total bandwidth of 250 MHz, cor- 
responds to 650 km s~^ at the ^^C0(1 - 0) frequency. 
The line intensity calibration was accomplished by the 
chopper- wheel method (Ulich & Haas 1976), and the rel- 
ative intensity errors were calibrated at every observa- 



tion with the intensity of W51 main for NGC 6951 and 
IRC+10216 for the rest of the observations. 

We observed the central regions of sample galaxies 
pointed toward the nuclei except for M82. M82 was ob- 
served toward five points for the ^^C0(1 - 0) and HCN(1 
- 0) lines, and three points for the ^^C0(1 - 0) line. The 
observed positions are shown in Fig. 2 and Table 3. The 
positions (0", 0"), (-15", -3"), and (10", 6") are the 
so-called the nucleus, the South- West (SW) lobe, and the 
North-East (NE) lobe, respectively (e.g., Shen & Lo 1995). 

The data were reduced using the reduction software 
package NEWSTAR of the Nobeyama Radio Observatory. 
The baselines of raw spectra were subtracted by fitting lin- 
ear lines, or if necessary, second order polynomials. The 
line-free velocity range for each galaxy was defined with 
avoiding the line emission velocity range based on the pre- 
viously published CO observations. The raw spectra with 
higher order baseline fluctuation were not used in our final 
spectra. 

4. Results 

The parameters of the observed lines from all of the 
sample galaxies (peak main beam brightness tempera- 
ture, Tmb, and integrated intensity, /) are summarized 
in Table 3, and the calculated integrated intensity ratios 
are summarized in Table 4. Since we did not observe the 
^^C0(1 - 0) line for the starburst galaxy samples other 
than M82 and the HCN(1 - 0) line for the non-starburst 
galaxy samples, we made use instead of the data taken 
with the Nobeyama 45 m telescope from various published 
papers (see notes in Table 4). 

Fig. 2 shows the "CO, "CO, and HCN line spectra 
toward the observed positions of the central region of 
M82, one of the starburst galaxy samples. The overall 
line shapes of each positions are similar across these three 
line spectra. The line shapes of all three spectra at the 
nucleus indicate constant brightness temperatures around 
the line center, Vi^sR ^ 150 — 300 km s^^. It is also obvi- 
ous that the ^^CO line is weaker than other two lines. On 
the other hand, the line shapes at the off-center indicate 
peak structures, and the peak brightness temperatures of 
the NE and SW lobes exceed those of the nucleus in all 
the observed lines. We did not detect the HCN line in the 
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^^CO(I-O) (0",0") 


A (-15",-3") 


(-30",-6") 









NE Lobe Nucleus SW Lobe 






0.8 


(10",6") 


HCN{1-0) (0",0") 


(-15" -3") 


(-30" -6") 


(-45" -9") 


Tmb [K] 


0.6 
0.4 
0.2 














0.0 

























200 400 200 400 200 400 200 400 200 400 
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Fig. 2. l2co(l - 0) (top row), "co(l - 0) (middle row), and HCN(1 - 0) (bottom row) spectra of the central region of the 
prototypical starburst galaxy M82. The horizontal axis is the LSR velocity, Vlsr, and the vertical axis is the main beam temperature, 
Tmb- The reference point (0", 0") is the position of the galactic nucleus as determined from the peak of the 2.2 /.tm source (see 
Table 1). The other points are in R.A. and Decl. offset from the reference point. 



farthest observed region, (—45", —9"). The line intensity 
ratios are, on the other hand, different from the brightness 
temperature distributions. At the center of M82, the in- 
tensity ratios (HCN/"GO, i^CO/^CO, and HCN/i^CO) 
are all high, but these ratios gradually decrease with in- 
creasing offsets from the center (Table 4). 

Fig. 3 shows 13C0(1 - 0) and HCN(1 - 0) line spectra 
of the other starburst galaxy samples. Both lines have 
similar line shapes and intensities for all galaxies, which 
is similar to those for the central ±15" region of M82. We 
note again that since the HGN and ^'^CO lines were ob- 
served simultaneously, there is no HGN/ -^"^ CO integrated 
intensity ratio errors due to the relative pointing offsets. 

Fig. 4 shows the "C0(1 - 0) and 1^00(1 - 0) line spec- 
tra of the non-starburst galaxy samples. The integrated 
intensity of ^^C0{1 - 0) line, l^^CO), for NGC 4736 is 
consistent with that of Kohno (1998). For 1{^^C0) of 
NGG 4826, it is 67% of the value obtained by Kohno 
(1998). This disagreement may be explained by point- 
ing error and/or calibration error. The pointing error of 
our observations was ^ 4", larger than the ^ 2" of Kohno 
(1998). On the other hand, Kohno (1998) did not observe 
any intensity calibrator as we did with IRG+10216. For 



I(i2C0) of NGC 5195, we only detect 46% of the value 
obtained by Kohno et al. (2002). This is because we had 
mistakenly placed the off-position of the position switch- 
ing observations at NGC 5194, which is the counterpart 
of the M51 interaction system for NGG 5195 and rich in 
molecular gas. This can be seen as the absorption-like 
feature around the velocity range of 500 — 550 km s~^ in 
Fig. 4, which did not show up in the spectrum of Kohno 
ct al. (2002). Therefore we use the derived HGN/"CO 
value for NGC 5195 as an upper limit. 

The HCN/^'^CO integrated intensity ratios between the 
starburst and non-starburst galaxy samples are signifi- 
cantly different: We obtained an average HCN/^"^GO to- 
ward the nuclei of the starburst samples of 1.15 ±0.32, 
and a value of 1.10 ±0.30 for starburst samples including 
the off-center starburst regions of M82, namely the NE 
and SW lobes. The average HGN/ "CO ratio toward the 
nuclei of the non-starburst galaxy samples and that in- 
cluding the non-starburst regions of M82, namely (—30", 
-6") and (-45", -9"), are < 0.31 ±0.14 and < 0.32 ±0.12, 
respectively. These values are consistent with previously 
published results: Paper I suggests that starburst regions 
tend to have higher HGN/ ^'^ CO of > 1, but less active star 
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Table 3. Observed properties for sample galaxies. 



Galaxy 




Position* 


TMB(^"CO)t 


TMB(^^CO)t 


rMB(HCN)t 






/(HCN)^ 


NGC 253 




(0", 0") 


— 


0.43 ±0.03 


0.42 ±0.02 


— 


168 ±15 


258 ±22 


NGC 2146 




(0", 0") 


— 


0.074 ±0.015 


0.072 ±0.007 


— 


16.0±1.3 


12.3±0.7 


NGC 2903 




(0", 0") 


— 


0.074 ±0.008 


0.066 ±0.009 


— 


9.1±0.7 


9.8 ±0.6 


M82 




(0", 0") 


3.52 ±0.08 


0.23 ±0.03 


0.27±0.02 


565 ± 33 


34.5 ±2.6 


55.4±2.6 


(NE lobe) 




f 10", +6") 


— 


0.76 ±0.04 


0.70 ±0.04 


— 


66.1 ±4.4 


71.2±3.6 


(SW lobe) 




-15", -3") 


6.02 ±0.13 


0.49 ±0.03 


0.39 ±0.03 


646 ± 38 


45.2 ±2.9 


36.2±1.9 






-30". -6") 


2.61 ±0.33 


0.27 ±0.05 


0.17±0.02 


288 ± 20 


29.0 ±2.8 


12.2±0.9 






-45", -9") 




0.17 ±0.05 


< 0.08 




12.7±2.2 


<3.2 


NGC 3504 




(0", 0") 




0.052 ±0.005 


0.060 ±0.009 




7.7±0.5 


10.1±0.6 


NGC 6946 




(0", 0") 




0.152 ±0.012 


0.151 ±0.008 




21.9 ±1.4 


21.4±1.0 


NGC 6951 




(0", 0") 




0.055 ±0.009 


0.034 ±0.007 




8.1±0.6 


6.1 ±0.5 


NGC 4736 




(0", 0") 


0.26 ±0.01 


0.043 ±0.012 




44.1 ±2.7 


6.3 ±1.4 




NGC 4826 




(0", 0") 


0.53 ±0.05 


0.11 ±0.03 




87.0 ±5.9 


17.2 ±3.6 




NGC 5195 




(0", 0") 


0.26 ±0.04 


0.10 ±0.02 




37.9 ±2.8 


12.0±2.2 





Notes. 

* Offset from the position of the nucleus indicated in Tables 1 and 2. 

t Peak main-beam temperature in K. Quoted uncertainties are la and upper limits are 2a. 

Integrated intensity in K km s~^. We define the integrated intensity as 7 = Jt^b'^'"- Quoted uncertainties are la and upper 
limits are 2a. 



Table 4. Integrated intensity ratios for sample galaxies. 



Galaxy 


Position* 


HCN/^^CQT 


^"CO/^^CQT 


HCN/^^CO* 


NGC 253 


(0", 0") 


1.53±0.19 


7.4 ± 0.611 


0.209 ±0.01811 


NGC 2146 


(0", 0") 


0.77 ±0.08 


11.8±1.0# 


0.065 ±0.004# 


NGC 2903 


(0", 0") 


1.08±0.11 


11.3±0.9** 


0.095 ±0.006** 


M82 


(0", 0") 


1.61±0.14 


16.4±1.6 


0.098 ±0.007 


(NE lobe) 


(±10", ±6") 


1.08 ±0.09 






(SW lobe) 


(-15", -3") 


0.80 ±0.07 


14.3±1.2 


0.056 ±0.004 




(-30", -6") 


0.42 ±0.05 


9.9±1.2 


0.042 ±0.004 




(-45", -9") 


< 0.25 






NGC 3504 


(0", 0") 


1.32±0.12 


8.8 ±0.6" 


0.094 ±0.008tt 


NGC 6946 


(0", 0") 


0.98 ±0.08 


9.7±0.6# 


0.101 ±0.005# 


NGC 6951 


(0", 0") 


0.75 ±0.08 


6.0 ±0.5** 


0.125 ±0.010** 


NGC 4736 


(0", 0") 


<0.32^ 


7.0±1.6 


< 0.042* 


NGC 4826 


(0", 0") 


0.48 ±0.10* 


5.0±1.1 


0.063* 


NGC 5195 


(0", 0") 


<0.13±0.06§ 


3.2 ±0.6 


0.018 ±0.007§ 



* Offset from the position of the nucleus indicated in Tables 1 and 2. 

t Ratio of integrated intensity indicated in Table 3. Integrated intensity that is not 
in Table 3 is taken from previously published papers, as indicated below. Quoted 
uncertainties are la, and upper limits are 2a, except for the HCN/^'^CO value of 
NGC 5195 (see note § below). 

t HCN(1 - 0) integrated intensity and HCN/l^CO ratio taken from Kohno (1998). 
§ HCN(1 - 0) integrated intensity and HCN/^^CO ratio taken from Kohno et al. 
(2002). Due to the off-position problem, we take HCN/^^CO as an upper limit (see 
Sect. 4 for detail). 

II ^^CO integrated intensity taken from Sorai et al. (2000). 

* -"^^CO integrated intensity taken from Sorai et al. (2002). 
** ^■^CO integrated intensity taken from Kuno et al. (2007). 

^^CO integrated intensity taken from Kuno et al. (2000). 
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Fig. 3. ^'^CO(l — 0) and HCN(1 - 0) spectra toward the center of the starburst galaxy samples. Two spectra are shown for each 
galaxy; the spectrum on the left-hand side denotes the ^''CO(l - 0) line, and that on the right-hand side denotes the HCN(1 - 0) 
line. The horizontal and vertical axes are Vlsr ^i^d Tmb i respectively. Note that the scale for the vertical axis is the same for the 
13 CO and HCN spectra. 



forming regions tend to have lower ratios (^1). Note that 
the results of Paper I is based on interferometric obser- 
vations; therefore it is sensitive to compact components, 
such as high density gas (i.e., HCN). It follows that ratios 
in Paper I are bound to take on larger values than those 
in this paper, which is based on single-dish observations. 

The i^cO/i^CO and HCN/i^CO integrated intensity 
ratios between starburst and non-starburst galaxy sam- 
ples also show differences with the average values for the 
starburst samples larger than those for the non-starburst 
samples: The average ^^CO/^^CO integrated intensity ra- 
tios toward the nuclei of the starburst and non-starburst 
samples are 10.2 ± 3.2 and 5.1 ± 1.6, respectively, and the 
average HCN/^^CO integrated intensity ratios for these 
samples are 0.11 ± 0.04 and < 0.04 ± 0.02, respectively. 

AU the HCN/"CO, i^cO/i^CO, and HCN/i^CO ratios 
for the starburst galaxy samples are 2 — 3 times larger 
than those for the non-starburst galaxy samples. These 
results suggest that the conditions of molecular gas in the 
starburst regions differ significantly from those in the non- 
starburst regions. 



5. Discussions 

5.1. Physical Conditions of Molecular Gas in Sample 
Galaxies 

Using the line ratios derived above, we can estimate 
the physical conditions of molecular gas in our starburst 
and non-starburst sample galaxies. To that end, we make 
use of the large-velocity-gradient (LVG) approximation 
(Goldreich & Kwan 1974; Scoville & Solomon 1974) as- 
suming a one-zone model, and the photon dominated re- 
gion (PDR) model (Meijerink & Spaans 2005; Meijerink et 
al. 2007). Note that the estimated physical conditions are 
averaged values within one beam size of our observations, 
which corresponds to a few hundred to a few kpc scale in 
linear size. Although deriving an averaged value within 
one beam is extremely simplified assumption compared to 
the real situations of molecular clouds and radiation fields 
in galaxies, it is successfully derived in various galaxies 
in the past studies. The overlapping effect of molecular 
clouds is expected to be small, based on the small shadow- 
ing effect in molecular clouds of our Galaxy that is edge-on 
(Burton & Gordon 1978), so we assume that we see most 
of the emission from each molecular clouds in our beam, 
and therefore it is safe to use the LVG and PDR models 
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Fig. 4. l^CO(l - 0) and ^^CO{l - 0) spectra toward the cen- 
ter of non-starburst galaxy samples. Two spectra are shown 
for each galaxy; the spectrum on the left-hand side denotes 
the ^^CO(l - 0) line, and that on the right-hand side denotes 
the ^^^00(1 - 0) line. The horizontal and vertical axes are 
Vlsr and Tmb, respectively. 

in our datascts. 

Spergel & Blitz (1992) suggested that the increase of the 
HCN intensity toward the centers of galaxies can be due to 
the increase of gas pressure, not due to the increase of gas 
density induced by gravitational instability. Hence, they 
predicted that larger galaxies should have stronger HCN 
intensities than those in smaller galaxies. This prediction 
can be rejected from the comparison of the HCN/^^CO 
and HCN/i^co ratios between our Galaxy and M82: Our 
Galaxy, which is a larger galaxy, has the HCN/ -'^^ CO and 
HCN/i^CO ratios of - 0.08 (Jackson et al. 1996; Heifer 
& Blitz 1997) and 0.8 (Matsushita et al. 1998), respec- 
tively, at the Galactic center. The center of a much smaller 
galaxy, M82, has much larger HCN/i^CO and HCN/i^CO 
ratios of ^ 0.1 and ~ 1.6, respectively. Similar conclusion 
is also obtained by Gao & Solomon (2004b) . We therefore 
assume that the differences of ratios between galaxies and 
within galaxies are due to the differences in physical con- 
ditions of molecular gas, not due to the differences in gas 
pressure. 

5.1.1. LVG Model 

The line intensity ratios for CO and HCN molecules 
under the LVG model were tabulated as a function of the 
kinetic temperatures, Tk, from 10 K to 1000 K, and the 



H2 number density, n(R2) 
for HCN(1 - 0)/i3CO(l - 0), 
and HCN(l-0)/i2CO(l-0). 



from 10^ cm-3 to 10^ cm'^ 
i2C0(l - 0)/i3CO(l - 0), 
The collision rates for CO 



molecules < 250 K and > 500 K were taken from Flower 
& Launay (1985) and McKee et al. (1982), respectively. 
The corresponding values for HCN molecules < 100 K and 
> 100 K arc available from Green & Thaddcus (1974) and 
on-line*. We fixed the abundances to the 'standard' 
relative abundance where Z{^^CO) = ["C0]/[H2] = 1 x 
10-^' (Solomon et al. 1979), Z(HCN) = [HCNj/pa] = 
2 X 10-8 (Irvine et al. 1987), and [i2CO]/[i3CO] = 50. 
We also set the velocity gradient, dv/dr, to a fixed general 
value of 1.0 km s"^ pc~^. The tabulated results arc shown 
in Fig. 5. 

The LVG calculations show that high HCN/i^CO ra- 
tio 1) can be obtained under high density (n(H2) ~ 
10^^^ cm"^) and high temperature (Tk ^ 50 K) condi- 
tions. The calculations for ^^CO/^'^CO show that high 
ratio (> 10) can be obtained under high density (n(H2) ~ 
10''='=^ cm"^) and high temperature (Tk > 100 K) con- 
ditions, or low density (n(H2) ^ 10^ cm^'^) conditions. 
The HCN/^^CO ratio, on the other hand, only traces the 
density, and a high ratio (^ 0.1) indicates high density 
(n(H2) > 2 X 10^ cm^'^) condition. These results can be 
explained as follows: At high temperature, the energy 
distribution of the "'^'^CO rotational transition populating 
higher levels due to the low abundance and small dipole 
moment of the ^^CO molecules, and therefore the ^^C0(1 
- 0) line will be optically thin with a corresponding de- 
crease in the line intensity. On the other hand, ^^CO 
molecules are abundant and HCN molecules possess large 
dipole moments, therefore the ^^C0(1 - 0) and HCN(1 - 
0) lines remain optically thick even under high tempera- 
ture conditions. In terms of density, the HCN molecule 
has a high critical density, meaning that high density is 
a prerequisite condition for line emission. The ^^CO and 
^•^CO molecules have the same critical density but ^-^CO 
emits under lower density conditions due to its abundance. 

^^CO deficiency is often cited in the case of merging 
galaxies (i.e., extreme starburst galaxies; e.g., Aalto et 
al. 1991; Casoli ct al. 1992). Our sample galaxies do not 
exhibit such extreme starburst behaviors, but the possi- 
bility of low Z{^^CO) cannot be ruled out. In addition, 
starburst galaxies may have larger velocity gradients due 
to the large turbulences caused by active star formation. 
In this study, we checked the dependence of Z(^^CO) and 
the velocity gradient with our LVG calculation results: If 
Z(^^CO) is a factor of two lower than the standard abun- 
dance, or if the velocity gradient is a factor of two higher 
than the standard value we assumed above, the LVG cal- 
culations for i^cO/i^CO and HCN/i^CO give rise to a 
temperature a factor of two lower. 

Using the LVG calculations shown in the figure, we 
arrived at the physical conditions of molecular gas in 
the sample galaxies and summarized them in Table 5. 
Many of the starburst galaxy samples display high den- 
sity (> 5 X 10"^ cm^-^) and high temperature (> 100 K). 
But some of the galaxies (NGC 2146 and NGC 6951) have 
been estimated as rather lower density or temperature due 



The data are available at 

http: / /data. giss.nasa.gov/mcratcs/data/hcn_he_ratcs.txt. 
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Fig. 5. H2 number density (vertical axis) and kinetic temperature (horizontal axis) dependence of the HCN(1 - 0)/^'^CO(1 - 0) 
(left), i2cO(l - 0)/13CO(1 - 0) (middle), and HCN(1 - 0)/12CO(1 - 0) (right) intensity ratios derived using the LVG calculations. 
Refer to the main text for the details of the calculations. 



Table 5. Physical conditions of molecular gas 



in the 


sample galaxies 


based LVG calculations. 


Galaxy 


Position* 


n(H2)T 




NGC 253 


(0", 0") 


- (1-2) X lO'' 


100 - 200 


NGC 2146 


(0", 0") 


- (3 - 5) X 10^ 


100 - 200 


NGC 2903 


(0", 0") 


- (6 - 9) X 10^ 


100 - 200 


M82 


(0", 0") 


- (6 - 9) X 10^ 


200 - 300 


(NE lobe) 


(+10", +6") 


-4x 104±i 


>50 


(SW lobe) 


(-15", -3") 


- (3-5) X 10^ 


-200 




(-30", -6") 


- (2-3) X 10^ 


- 100 




(-45", -9") 


< 5 X 10^ 




NGC 3504 


(0", 0") 


- 1 X 10* 


-200 


NGC 6946 


(0", 0") 


- (6 - 9) X 10^ 


100 - 200 


NGC 6951 


(0", 0") 


- (6 - 9) X 10^ 


60 - 100 


NGC 4736 


(0". 0") 


< 2 X 10^ 


^60 


NGC 4826 


(0", 0") 


- (3 - 4) X 10^ 


- 50 - 60 


NGC 5195 


(0". 0") 


- 1 X 10^ 


- 10-20 



INotes. 

* Offset from the position of the nucleus indicated in Tables 1 and 2. 
t Molecular hydrogen number density in cm~^. 
Kinetic temperature in K. 



to lower HCN/^^CO values, even these galaxies have high 
infrared luminosities (see Table 1). This is probably be- 
cause of the larger distances of these galaxies than other 
nearby galaxies: Due to larger distances, the beam ob- 
serves larger radius, where usually dominated by low den- 
sity and low temperature materials, than other nearby 
galaxies, and therefore the beam averaged physical condi- 
tions exhibit lower values than other nearby galaxies. On 
the other hand, the non-starburst galaxy samples displays 
lower densities and temperatures of n(H2) ^ 4 x 10'^ cm~'^ 
and Tk ^ 60 K than the starburst galaxy samples, even 
non-starburst samples located nearby. 

We then compared our estimations with previously pub- 
lished values. In the case of M82, a region of higher density 
(—4 X 10"'^"'^ cm~'^) can be observed at the NE lobe, which 
gradually tapers off toward the southwestern side of the 



galaxy. This trend, as well as the corresponding derived 
values, are very similar to those derived by Petitpas & 
Wilson (2000). The temperature, on the other hand, has 
high values (200 - 300 K) at the nucleus and the SW lobe, 
and — 100 K at the position (—30", —6"). The decrease 
of the temperature as a function of radius in M82 (and 
in other galaxies) is also suggested from the ^^CO/^'^CO 
observations (Paglione et al. 2001). The temperature of 
the NW lobe is not well constrained, ^ 50 K, and there 
exists no constraint for the position (—45", —9") due to 
the absence of ^^C0(1 - 0) data. These estimated temper- 
atures and densities tend to be higher than those derived 
using various CO lines (including isotopes) observed with 
an interferometer (Weifi et al. 2001). 

Estimations for NGC 253 yield n(H2) - (1 - 2) x 
10* cm~'^ and Tk — 100 — 200 K, which are very similar 
to previously published results, utilizing several molecular 
species (Wall et al. 1991; Paglione et al. 1995; Bradford 
et al. 2003; Bayet et al. 2004), but rather higher than 
those derived using only *^C0 lines (Glisten et al. 2006). 
For NGC 6946, our estimations suggest that n(H2) — 
(6 - 9) x 10^ cm-3 and Tk - 100 - 200 K, which are simi- 
lar to the values obtained by Meier & Turner (2004) using 
various CO lines including isotopes and also HCN(l-O) 
line, but rather higher than those derived by Walsh et al. 
(2002), which only made use of the CO lines (including 
isotopes) . 

The above comparisons indicate that our approach 
tends to yield higher densities and temperatures than 
other studies using only the ^^CO lines or even the CO iso- 
topes. This is because our methodology makes extensive 
use of both the high density tracer HCN(l-O) line and 
the optically thin and therefore temperature dependent 
tracer 13CO(1-0) line. In particular, HCN(l-0)/i3CO(l- 
0) is sensitive to both density and temperature, once it 
attains a value of around unity or higher (see Paper I 
for more detailed discussions). We therefore concluded 
that the combinations of these three line ratios, or even 
HCN(l-0)/i3cO(l-0) by itself, can effectively distinguish 
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Table 6. IRAS data and infrared flux ratios of sample galaxies. 



[Vol. , 



. J(12 Mm) 
/(25 Mm) 



/(12 Mm) 
/(60 Mm) 



Mm) 
/(lOO Mm) 



^ /(25 Mm) 
^" 6 /(60 Mm) 



"7(25 



/(60 Mm) 
/(lOO Mm) 



Galaxy 



lo! 



log 



log 



/(lOO Mm) 



log 



NGC 253 
NGC 2146 
NGC 2903 

M82 
NGC 3504 
NGC 6946 
NGC 6951 



-0.5762 ± 
-0.4400 ± 
-0.2131 ± 
-0.6220 ± 
-0.5600 ± 
-0.2328 ± 
-0.2073 ± 



0.0004 
0.0016 
0.0045 
0.0002 
0.0125 
0.0013 
0.0053 



-1.3726 ±0.0001 
-1.3320 ±0.0005 
-1.0586 ±0.0018 
-1.2704 ±0.00005 
-1.2857 ±0.0034 
-1.0301 ±0.0007 
-1.0835 ±0.0029 



-1.4968 ±0.0004 
-1.4535 ±0.0015 
-1.3919 ±0.0029 
-1.2379 ±0.0003 
-1.4868 ±0.0116 
-1.3803 ±0.0013 
-1.4938 ±0.0049 



-0.7964 ±0.0001 
-0.8920 ±0.0007 
-0.8455 ±0.0035 
-0.6484 ±0.00007 
-0.7257 ±0.0049 
-0.7972 ±0.0007 
-0.8761 ±0.0029 



-0.9206 ±0.0003 
-1.0135 ±0.0009 
-1.1789 ±0.0035 
-0.6159 ±0.0002 
-0.9268 ±0.0049 
-1.1475 ±0.0009 
-1.2864 ±0.0034 



-0.1242 ±0.0002 
-0.1215 ±0.0006 
-0.3333 ±0.0008 
+0.0325 ±0.0002 
-0.2011 ±0.0015 
-0.3502 ±0.0007 
-0.4103 ±0.0020 



NGC 4736 
NGC 4826 
NGC 5195t 



-0.0810 ±0.0055 
-0.0835 ±0.0105 
-0.4645 : 



-1.1495 ±0.0040 
-1.1918 ±0.0048 
-1.6383 : 



-1.3767 ±0.0048 
-1.5390 ±0.0057 
-1.9512 : 



-1.0685 ±0.0029 
-1.1083 ±0.0089 
-1.1738 : 



-1.2956 ±0.0029 
-1.4556 ±0.0089 
-1.4874 : 



-0.2271 ±0.0009 
-0.3473 ±0.0011 
-0.3135 : ±0.0228 



Notes. 

Flux data are taken from Sanders et al. (2003). 

t A value with a colon (":") indicates that the value has a large uncertainty. 



galaxies with warm and dense molecular gas from galaxies 
without. 

5.1.2. PDR Model 

We also calculated HCN(l-0)/i3CO(l-0) with the PDR 
model. It turns out that the PDR model can also yield 
a HCN/^'^CO of around unity; it requires an interstellar 
radiation field (ISRF) strength of around several xlO'^ Gq 
and a density of around several xlO^ cni""^. In other 
words, dense molecular gas irradiated by radiation from 
massive stars can reproduce a high HCN/^'^CO of around 
unity. The ratio decreases with decreasing ISRF strength 
or density, but the decrease of the ratio is sensitive to the 
density decrease than to the ISRF strength decrease. 

Using this PDR model, the ISRF strength and the den- 
sity for the starburst sample are estimated as around 
10^"'^ Go and on the order of 10^ cm"-^, respectively. The 
non-starburst sample, on the other hand, has lower ISRF 
of the order of 10^ Go and lower density of the order of 
10'* cm^'^ or less. 

The strength of the ISRF can be derived also from 
the IR characteristics. The range of the ISRF strength 
in nearby galaxies is ^ 10^~^ Go (Negishi et al. 2001; 
Malhotra et al. 2001), and exhibits a tight correlation 
with the IR flux ratios between 60 fim. and 100 /im, 
/(60 /.tm)//(100 /.tm); higher the ratio, higher the ISRF 
strength (Negishi et al. 2001). Our sample galaxies (both 
starburst and non-starburst samples) have a wide range 
of /(60 /im)//(100 A*in)i f^nd our starburst sample has a 
tendency of higher ratio than that of the non-starburst 
sample, although the range is overlapping with the non- 
starburst galaxy sample (see Fig. 1 and Table 6). If we 
adopt the correlation shown by Negishi et al. (2001), our 
starburst sample has higher ISRF strength of the order 
of 10^"'^ Go than the non-starburst sample of the order 
of 10^ Go. Those values are consistent with the values 
derived from HCN/^^CO using the PDR model, suggest- 
ing that there is a correlation between HCN/^^CO and 
/(60 /zm)//(100 /im). In the next subsection, we study 
this relation, and the relations between HCN/^^CO and 
other IR flux ratios. 



5.2. Relation between Molecular Gas and Dust 

Dust in galaxies absorbs the radiation from stars, in 
particular massive stars, and re-radiate at IR wavelengths. 
The IR emission from galaxies, especially IR flux ratios, is 
therefore a good tracer for star formation activities (e.g., 
Helou 1986; see also reviews of Soifer et al. 1987 and 
Telesco 1993). The spectral energy distribution study us- 
ing the IRAS and the Infrared Space Observatory (ISO) 
data between 3 and 1100 /xm further confirmed the utility 
of the IR flux ratios for tracing star formation activities, 
and also suggested that the mid-IR continuum flux at a 
range of 20 - 42 /xm can be the optimum dust emission 
tracer for current star formation activities (Dale et al. 
2001). Here we compare HCN/^'^CO and IR continuum 
flux ratios to ascertain the relation between molecular gas 
and dust properties. 

In this study, we use the IRAS flux (Sanders et al. 
2003) for the infrared flux due to the following reasons: 
(1) The IRAS data are still the only all-sky survey results 
so far, and all the wavelength data exist for all our sam- 
ple galaxies. (2) Although the beam size of the IRAS 
data is large (arcminutes scale), most far-IR radiation 
from nearby galaxies is concentrated in the central ^ 30" 
(Smith & Harvey 1996). The regions that traced by our 
observations and that by the IRAS are therefore more 
or less similar. Indeed, previously published single dish 
molecular gas studies derived important results using the 
IRAS data. The IRAS flux ratios for our sample galaxies 
are tabulated in Table 6. 

Since /(12 /im)//(25 /im) and /(60 /im)//(100 /im) 
arc good tracers of star formation (see Sect. 2 and 
Fig. 1), and the relation between HCN/^^CO and 
/(60 /im)//(100 /im) is suggested in the previous sub- 
section, we first compare HCN/^^CO with these IR ratios 
(Fig. 6a, b). HCN/^^CO indeed exhibits a correlation 
with /(60 /im)//(100 /tm). There exists also a correla- 
tion between HCN/^^co and /(12 /im)//(25 /im) ratio, 
which is natural since /(12 /im)//(25 /im) correlates with 
/(60 /im)//(100 /im) (Helou 1986; Dale et al. 2001). In 
addition, the correlation is much tighter than the case 
with /(60 /im)//(100 /im), except for NGC 5195 that has 
a large uncertainty in the IR flux (see Table 6). 
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Fig. 6. Correlation diagrams between the HCN{1 - 0)/^''CO(1 - 0) integrated intensity ratios and (a) /(60 ^m)//(100 /jm), (b) 
/(12 fim)/ f{25 fim), (c) /(25 /im)//(60 fim), and (d) /(25 /xm)//(100 fim) IR flux ratios. FiUed and open circles are the same as 
in Fig. 1. Downward arrows denote upper limits. All the diagrams, especially (b), (c), and (d), exhibit tight correlations. 



In Fig. 6(c) and (d), we also included correlation dia- 
grams for /(25 /xm)//(60 fj,m) and /(25 /im)//(100 /^m), 
and both show tight correlations as seen in the HCN/^'^CO 
- /(12 jim)/ f{25 fim) diagram. On the other hand, there 
is no significant correlation with other IR flux ratios, 
namely /(12 /xm)//(60 fim) and /(12 fim)/ f {100 fim) (di- 
agrams not shown). 

The tight correlations between HCN/^'^CO and the 
25 fim related IR flux ratios indicate that the increase 
of the 25 fim continuum flux relative to the other IR con- 
tinuum bands is tightly correlated with the increase of 
HCN/"CO. Under the LVG model, HCN/i^CO is sensi- 
tive to molecular gas that has a temperature of 100 K or 
more. Under the PDR model, HCN/^^cq jg sensitive to 
the molecular gas density and the strength of the ISRF. 
The increase of the 25 fim continuum flux is mainly due 
to the increase of flux from very small grains (Dale ct 
al. 2001). It is believed that the heating of these grains 
is intermediate between thermal equilibrium and single- 
photon heating, and the increase in flux is closely related 
to the increase in the heating intensity of the ISRF (Dale 
et al. 2001). This is consistent with the results from the 



PDR model. Furthermore, the black body temperature 
around 25 fim is around 100 K or more, therefore dust 
needs to be heated up to approximately this temperature 
before achieving thermal equilibrium, in complete agree- 
ment with the results from the LVG model. The tight 
correlations therefore suggest that the molecular gas and 
dust are localized in the same regions, with the strong 
radiation from starburst activities heating the dust and 
in turn the molecular gas. With the dust heated up to 
around 100 K or more, the molecular gas will follow suit. 
In cither LVG or PDR model, this tight correlation sug- 
gests that HCN/^'^CO will be an excellent tracer for star- 
burst conditions. It also applies that either or both the 
LVG and PDR mechanisms play important roles in star 
forming regions. 

It is now well established that the amount of dense 
molecular gas is closely correlated to star formation ac- 
tivities (Solomon et al. 1992; Gao & Solomon 2004a; Gao 
& Solomon 2004b). Our results suggests that in addi- 
tion to the amount of dense molecular gas, the amount of 
dense and warm molecular gas is also pertinent to star for- 
mation activities. The possible causes of this correlation 
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are the heating of molecular gas due to strong radiation 
and stellar wind from newly formed massive stars, and/or 
numerous supernova explosions. 

6. Conclusion 

We conducted a survey of ^^C0{1 - 0), ^^00(1 - 0), 
and HCN(1 - 0) lines toward seven starburst and three 
non-starburst galaxies nearby using the Nobeyema 45 m 
telescope (beam size ~ 14" - 18"). The ^^Coll - 0) and 
HCN(1 - 0) lines were obtained for all the sample galaxies, 
but the ^^C0(1 - 0) data for starburst galaxy samples 
and the HCN data for non-starburst galaxy samples were 
obtained from previously published data. 

• The i3C0(l - 0) and HCN(1 - 0) intensities in the 
central regions of the starburst galaxy samples are 
more or less similar, with the HCN/^^CO integrated 
intensity ratios ranging between 0.75 - 1.61 (average 
= 1.15 ±0.32). The ^^C0{1 - 0) intensities in the 
non-starburst galaxy samples are always stronger 
than the HCN(1 - 0) intensities, and HCN/i^CO 
yield obviously smaller values of < 0.31 ± 0.14 than 
for the starburst galaxy samples. We observed mul- 
tiple positions for the starburst galaxy M82, and the 
ratios far from the center (i.e., disk region) are sim- 
ilar to those in the non-starburst galaxy samples. 

• i^cO/i^CO and HCN/i^CO also show higher values 
for about a factor of two in starburst galaxy samples 
than those in the non-starburst galaxy sample. 

• The physical attributes of molecular gas derived 
from the line ratios under the LVG model indicate 
dense (> 5 x 10'^ cm~'^) and warm (> 100 K) con- 
ditions in most of the central regions of the star- 
burst galaxy samples, but diffuse (^ 4 x 10'^ cm"'^) 
and cold (^ 60 K) conditions in the non-starburst 
galaxy samples. Under the PDR model, this sug- 
gests that the dense molecular gas (^ 10^ cm~'^) is 
being irradiated by ISRF of 10^^"^ Go in the star- 
burst galaxy samples, but lower density molecular 
gas of ^ lO'' cm~^ is being irradiated by lower in- 
tensity ISRF of 10^ Go in the non-starburst galaxy 
samples. 

• HCN/ CO exhibits tight correlations with the 
IRAS flux ratios between 25 fim and other wave 
bands, namely, the increase of HCN/^'^CO is re- 
lated to the increase of the IRAS 25 fim flux. Since 
the IRAS 25 /.m flux is a good tracer for star for- 
mation activities, these correlations indicate that 
HCN/ ^"^ CO is also a good tracer for star formation 
activities. High HCN/^'^CO larger than unity points 
to dense and warm molecular gas with the tempera- 
ture of 100 K or more under the LVG model, or irra- 
diated by strong ISRF under the PDR model. The 
high 25 ^m flux also indicates high ISRF and high 
dust temperature of around 100 K or more. These 
results suggest that molecular gas and dust are lo- 
calized in the same regions in a well mixed configura- 
tion, and that the strong radiation heats up the dust 



and therefore the molecular gas nearby indirectly. 
This strong radiation also affects the molecular gas 
itself. 
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Taiwan, NSC 96-2112-M-001-009 and NSC 97-2112-M- 
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